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We give detailed predictions for the diffractive associated production of SUSY Higgs bosons and 
top squarks at the LHC via exclusive double pomeron exchange mechanism. We study how the SUSY 
Higgs production cross-section and the signal-over-background ratio are enhanced as a function of 
tan/3 in different regimes. The prospects are particularly promising in the "anti-decoupling" regime, 
\^~) ■ which we study in detail. We also give prospects for a precise measurement of the top squark mass 

' using the threshold scan of central diffractive associated top squark events at the LHC. 
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I. INTRODUCTION 

<\ 

The quest for supersymmetry (SUSY) is one of the major goals in high energy particle physics, and is an important 
task for the experiments at the upcoming LHC. The discovery of SUSY Higgs scalar (s) and other new SUSY particles, 
among which the superpartners of the top quark are good candidates, would be an important achievement. Standard 
production mechanisms based on QCD are now well explored, at least for the main channels. However, due to both 
the general interest of the problem and some specific features of the SUSY Higgs and top squark sectors, one should 
seek for alternative ways of SUSY production. In this paper we investigate the prospects for diffractive production of 
\Qi SUSY Higgs bosons and associated sparticles (stops) in the central region of the detectors. 

Standard Model Higgs boson production in double diffraction (denoted DPE, for Double Pomeron Exchange) has 
already been studied in recent years Q, 0, QjGLIlEl EE El Many approaches have been pursued, considering 
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diffractive scattering, as in the Regge picture PHSQi final state soft colour interactions 0, or fully perturbative 
exchange of gluon pairs [7j ■ We extend these studies to the SUSY Higgs and sparticle sector in the framework of the 
Minimal Supersymmetric Standard Model (MSSM) This subject has already been investigated in the literature 
JL • [l2l IT^j l 1 , and we will in particular focus on central diffractive production of the lightest MSSM Higgs boson in the 
' "anti-decoupling" regime which has not yet been studied in the proposed framework. The "decoupling" and "intense 
coupling" regimes have been studied in Ref. 0] > where it is shown that diffractive Higgs boson production can help 
L" ' distinguishing between h and H in the intense coupling regime. 
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S-h ' A. Central diffractive production of a heavy state 

One generally considers two types of DPE events for the production of a heavy state, namely "exclusive" DPE 
[JHHEl' where the central heavy object is produced alone, separated from the outgoing hadrons by rapidity gaps : 

pp — > p + heavy object + p , (1) 
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and "inclusive" DPE 0, 0, IE S 0i where the colliding Pomerons are resolved (very much like ordinary hadrons), 
accompanying the central object with Pomeron "remnants" (X,Y): 

pp — > p + X + heavy object + Y + p . (2) 

In general, exclusive production is considered most promising, because of a good expected signal-over-background 
ratio (due to the large gaps with no or low hadronic activity specific of diffractive events) and because of the good 
missing mass resolution [l4j 2 . Obviously, hard diffractive cross-sections are of higher order than standard hard 
non-diffractive ones, and this implies lower cross-sections. 

There are two objectives for this paper: First, to present a detailed calculation of the cross-section for exclusive 
DPE events, and second, to elaborate on the advantage of using these events over other processes to search for new 
heavy objects and to determine their characteristics. Indeed, if the events are exclusive, i.e., if no other particles 
are produced in addition to the pair of heavy objects and the outgoing protons, the measurement of the scattered 
protons in roman pot detectors allows to access the Pomeron-Pomeron centre-of-mass [14| , and to study accurately the 
dynamics of the hard process. It is therefore possible to measure with great accuracy the properties of new particles, 
e.g. their mass. It is also possible to study the new couplings by measuring cross-sections and angular distributions. 
As an example of this approach, we give a detailed description of the Higgs boson production cross-section and the 
top squark mass measurement at production threshold. The method can easily be extended to other heavy objects. 

These studies rely on the existence of exclusive events. At present, collider data from the Fermilab Tevatron 
give no evidence for the existence of exclusive production processes. Only upper bounds are given [15j, dominated 
by the study of dijet events near the kinematic limit of the produced diffractive mass. It is particularly difficult to 
identify purely exclusive events experimentally, since the pomeron remnants are not detected, and the mass resolution 
measured in the D0 or CDF main detectors suffer large uncertainties due to the detector resolution. The evidence 
for exclusive events would be much clearer in the 7 channel since the measurement of the diphoton mass fraction does 
not suffer from these uncertainties. However, the production cross section of those events is expected to be small at 
the Tevatron and LHC data will be necessary to observe this process. Testing the present model with Tevatron data 
will thus be challenging and LHC data in the beginning of data taking will be of fundamental importance. 

The theory of exclusive events is uncertain as well. Different models lead to cross-sections at the LHC differing 
by orders of magnitude. Two approaches are commonly discussed in the literature |l6|. hereafter referred to as the 
'Pomeron- induced' 0,0] and 'proton- induced' Q models 3 for exclusive production. They share the common feature of 
satisfying the upper bounds for exclusive production at the Tevatron and they give similar predictions for a low-mass 
standard Higgs boson cross-sections at the Tevatron. Apart from this, they come from different dynamics, since the 
'proton-induced' model is based on a semi-hard perturbative gluon mechanism, while the 'Pomeron-induced' model 
is based on a soft, essentially non-perturbative, mechanism based on Pomeron exchange. 

The 'Pomeron-induced' model is an extension Q to the purely exclusive processes of the original Bialas- 
Landshoff process for diffractive Higgs production which was applied to heavy quark pairs in 2]. In our model, 
both inclusive and exclusive diagrams come from the same approach and are beased on the Bialas Landshoff model 
The model starts with the same soft Pomeron exchange diagrams (for ordinary dijet production the 00 — > gg 
diagrams are also included 0) and correct the result by non-perturbative rapidity-gap survival factors The 
energy dependence is related to the rise of ordinary hadronic cross-sections through features of the soft Pomeron [2(| . 

An important issue of the 'Pomeron-induced' model is that the colourless Pomeron exchange implies the diffractive 
phenomenon. Our model follows the Bialas Landshoff model presciption, namely we assume the existence of a direct 
coupling between the Pomeron and the heavy mass object which is produced (Higgs or jets for instance). Hence, by 
definition it takes into account the veto on gluon radiation from the production mechanism. In other words, since we 
consider a pomeron induced model, there is no further need for a Sudhakov form factor to suppress radiation at the 
pomeron level. The counterpart is that it involves non perturbative mechanisms which have to be modelled using 
non perturbative gluon propagators. Hence, in the 'Pomeron-induced' framework, the normalisation is not determined 
theoretically, since it is related to the unknown non-perturbative strong coupling constant G 2 /4tt 0,0- However, the 
coupling constant value can be fixed phenomenologically by requiring consistency with the description of the inclusive 
dijet production within the same scheme, and is found Q compatible with the value G 2 /An — I chosen in 1, 2}. For 
our study, we will use this value to evaluate the production cross-sections 



2 The missing mass can be computed very precisely using roman pot detectors, and is equal to twice the mass of the heavy object in the 
case of exclusive events. 

3 We call 'proton-induced' the model of Ref.Q], since it has the unintegrated gluon distribution inside the proton as the main ingredient. 
This is contrasted with the 'Pomeron-induced' model 0, where the main source are gluons in the Pomeron. The 'proton-induced' 
and 'Pomeron-induced' models are sometimes called respectively, the 'Durham model' and the 'Saclay model' for exclusive production 
in the literature III! . 



The 'proton-induced' model is based on the perturbative calculation related to the suppression of QCD radiation 
due to Sudakov form- factors. Besides, soft radiation is forbidden through a rapidity gap survival formalism [Ls| . The 
same Sudakov form-factors are responsible for a damping of large mass diffraction (for Higgs boson production, it is 
compensated by a rapid growth as a function of energy). In particular, one expects to get a negligible top squark pair 
exclusive cross-section [21]. Indeed, in the 'proton-induced' model approach, the normalisation is estimated purely 
theoretically (except the lower bound constraints on dijet production) and leads to a top squark cross-section too 
small to be observed at the LHC. 

We will focus on exclusive production and restrict ourselves to the original Bialas-Landshoff type of models 
0,0,13, with their extension to SUSY Higgs and stop production which we will develop in the next sections. There 
are large theoretical uncertainties of diffractive production which are related to the non trivial interplay between 
perturbative and non perturbative contributions. We found it instructive to use the non perturbative 'pomeron 
induced' model and compare it to the predictions for the 'proton-induced' model. Indeed, it will be straightforward to 
extend these studies to other models. A recently developed Monte-Carlo program, DPEMC [22^. implements the models 
of [aaaiH. Moreover, most of the plots are in terms of s/b and enhancement factors, that are independent of 
the models. 

Even though inclusive DPE is a less promising search channel, it is still important to consider. In particular, 
it is interesting to evaluate the tail of the inclusive mass spectrum ("quasi-exclusive" processes) since it constitutes a 
background to exclusive DPE. In addition, only inclusive DPE has actually been observed for high central masses [lflj . 
Issues of inclusive DPE will be discussed in a future publication. 



B. The relevant SUSY spectrum: SUSY Higgs boson and top squarks 

Due to the limitation in the available total energy for production, diffractive production is favoured for the pro- 
duction of SUSY particles in the lower range of their mass in the admissible set of model parameters. Hence, we 
will focus on this range for the SUSY Higgs and top squarks sector. The regions of the MSSM parameter space that 
favour a light Higgs boson and those that favour a light top squarks are not the same, so a specific study is required 
separately for Higgs bosons and top squarks, see Section III and IV. 

It is well-known that the Higgs boson sector of the MSSM is richer than that of the Standard Model. First, it 
contains five physical scalar degrees of freedom, instead of a single one: two CP-even neutral Higgs bosons h and H, 
a pseudoscalar Higgs bosons A and a charged Higgs boson pair H . Secondly, the lightest MSSM Higgs boson h may 
look very different from the SM Higgs boson. One can define (at least) three noteworthy regimes for the couplings of 
the neutral MSSM Higgs bosons h, H and A: 

(i) the decoupling regime, in which h behaves like the SM Higgs boson [2^ll23|: 

(ii) the intense coupling re gim e, in which the couplings of all three neutral Higgs bosons are very different from 
those of the SM Higgs boson [2^] ; 

(iii) the so-called anti- decoupling regime |2fij . in which H behaves like the SM Higgs boson, while h has enhanced 
(resp. suppressed) couplings to down-type fermions (resp. up- type fermions and gauge bosons) |27j . 

It is also well-known that the MSSM Higgs boson sector contains at least one scalar h with rather low mass (the 
other, H , being with larger but possibly accessible mass) which gives a particular interest for diffractive production 
as we will study in detail in this paper. Indeed, the small mass, the sometimes small rate and the experimentally 
difficult standard channels, e.g., the decay into 77, enhances the interest in alternative production modes and decays 
such as diffractive production. We will in particular focus on central diffractive production of the lightest MSSM 
Higgs boson in the anti-decoupling regime (iii) which has not yet been intensively studied in the proposed framework. 
The decoupling and intense coupling regimes have been studied in Ref. [ill, where it is shown that diffractive Higgs 
boson production can help to distinguish between h and H in the intense coupling regime. 

The interest of MSSM Higgs boson production via exclusive diffractive production parallels a similar analysis for 
the Standard Model Higgs boson, with some distinctive features which enhance the specific production and branching 
modes. 

The discovery of "sparticles" at the LHC would be the clearest and most exciting signal of new fundamental physics 
beyond the Standard Model. Among these, the scalar superpartners of the top quark are expected to be those with 
smallest mass among scalars in a large portion of the MSSM parameter space. Indeed, various supersymmetric 
scenarios can accommodate a light top squark consistent with the experimental bounds on other sparticle masses 
and with measurements of observables that could be affected by large supersymmetric contributions, such as the 
anomalous magnetic moment of the muon or the branching ratio of the flavour violating decay b — > sj. Minimal 
supergravity [28j (mSUGRA) scenarios with a light top squark typically require a low gaugino mass parameter m 1 / 2 
and a large A-term. parameter Aq. The need for a small m±/2 is due to the fact that the renormalisation group 



equations for the soft supersymmetry breaking squark masses Mq and Mj^ receive a large contribution from gluinos: 
the larger m 1 / 2 , the higher the weak-scale values of Mq 3 and M^ 3 . As an example, the Snowmass Point 5 (SPS 5), 
defined by the following values of the mSUGRA parameters: m = 150 GeV, m 1/9 = 300 GeV, A n = —1000 GeV, 
tan/3 = 5 and sign(/i) = +, yields the following top squark and sbottom spectrum |29| . as generated by the program 
SUSYGEN 3.00/27 

m~ ti = 210 GeV , m t ~ 2 = 632 GeV , m gi = 561 GeV , mg 2 = 654 GeV . (3) 

For comparison, a "typical" mSUGRA scenario with vanishing Ao, the "post-WMAP benchmark scenario" B', defined 
by m = 60 GeV, m 1/2 = 250 GeV, A a = 0, tan/3 = 10 and sign(» = +, yields jp: 

m~ ti = 393 GeV , m h = 573 GeV , m gi = 502 GeV , mg 2 = 528 GeV . (4) 

Light top squarks and bottom squarks 4 can also arise from non-minimal SUGRA models, e.g. from scenarios with an 
inverted mass hierarchy in the squark sector [3^]. Assuming M| 3 <C a ($ = Q, J7, £)) at the GUT scale and small 
gaugino masses, one ends up with very low third generation squark masses at the weak scale due to strong 2-loop 
renormalisation group effects proportional to a|Tr(2Mg + + M|>) On the contrary, the first two squark 

generations remain heavy. The top squark and sbottom squared masses can even be driven negative if the GUT-scale 
hierarchy Mj 3 <C a is too pronounced, or if the gluino mass, whose contribution to the running of the squark 
masses tends to compensate for the two-loop gauge contribution, is too small. 

The paper is organised in the following way. In the next Section II, we introduce the concept of central diffractive 
production of SUSY Higgs bosons and top squarks; in II-A, the formalism of exclusive production and in II-B the 
experimental context are presented. In section III, we focus on the MSSM Higgs boson sector; in III-A, theoretical 
aspects of the Higgs boson spectrum and in III-B, the predictions for the LHC are displayed. In section IV, the 
case for top squark, and eventually bottom squark production is discussed; in IV-A, the theoretical framework, in 
IV-B, the predicted cross-sections and missing mass distribution and in IV-C, the top squark mass measurement by 
a threshold scan are given. The paper ends by a conclusion and outlook. 



II. DIFFRACTIVE PRODUCTION OF SUSY HIGGS BOSON AND TOP SQUARKS 

A. Exclusive central diffractive production 

In this section we introduce the model 0, 0, Q that is being used to describe exclusive MSSM Higgs bosons and top 
squark pair production in double diffractive production. In 0, , the diffractive mechanism is based on two-gluon 
exchange between the two incoming protons. The soft pomeron is seen as a pair of gluons coupled non-perturbatively 
to the proton. One of the gluons is then coupled perturbatively to the hard process, either the SUSY Higgs bosons, 
or the tt pair, while the other one plays the role of a soft screening of colour, allowing for diffraction to occur. The 
corresponding cross-sections for Higgs bosons and tt production read: 

dalf c (s) = <% 5^ I (vt + h) I II dm if vl exp(-2A tt ^ 2 )} a(gg - tt ) (5) 

^ tt 2 ' \i=l,2 J i=l,2 

where, in both equations, the variables vi and £j denote the transverse momenta and fractional momentum losses 
of the outgoing protons. In the second equation, fej and m are the squark transverse momenta and rapidities. 
a (99 ~ * H),a{gg — > tt ) are the hard production cross-sections which are given later on. The model normalisation 
constants Ch, C^ are fixed from the fit to dijet diffractive production, and are given in Ref . P, 0, H3| • 

In the model, the soft pomeron trajectory is taken from the standard Donnachie-Landshoff parametrisation |20| . 
namely a(t) = 1 + e + a't, with e 0.08 and a' ps 0.25GeV~ 2 . \h, Ay? are kept as in the original paper PJQ for the 



4 In the following, we consider only top squarks. But the study remains unchanged for squarks (i.e., bottom squarks) if their masses are 
sufficiently low. 



SM Higgs boson and qq pairs 5 . Note again that, in this model, the strong (non perturbative) coupling constant is 
fixed to a reference value G 2 /An = 1. 

In order to select exclusive diffractive states, one needs to take into account the corrections from soft hadronic 
scattering. Indeed, the soft scattering between incident particles tends to mask the genuine hard diffractive interactions 
at hadronic colliders. The formulation of this correction 0,^1 leads to an overall correction factor of 3% which is. 
the commonly used correction factor 0,0] for the QCD exclusive diffractive processes at the LHC. 

B. Experimental context 

The DPEMC [2^ Monte Carlo program provides an implementation of the Higgs boson, top squark and bottom 
squark pair production described above in both exclusive and inclusive double pomeron exchange modes. It uses 
HERWIG [34| as a cross-section library of hard QCD processes and, when required, convolutes them with the relevant 
pomeron fluxes and parton densities. The survival probabilities discussed in the previous section (0.03 for double 
pomeron exchange processes) have been introduced at the generator level. The cross-sections at the generator level 
are given in the next section after this effect is taken into account. 

A possible experimental setup for forward proton detection is described in detail in 0iH^|. We will only describe its 
main features here and discuss its relevance for the Higgs boson and top squark or bottom squark mass measurements. 

In exclusive DPE or QED processes at the LHC, the mass of the central heavy object can be reconstructed using 
the roman pot detectors and tagging both protons in the final state. It is given by M 2 — £i£2S, where £j are the 
proton fractional momentum losses, and s is the total centre-of-mass energy squared. 

In the following, we assume the existence of two detector stations, located at approximately 210 m and 420 m 
from the interaction point. The £ acceptance and resolution have been derived for each device using a complete 
simulation of the LHC beam parameters 35]. The combined £ acceptance is close to ~ 60% at low masses of about 
100 GeV, and 90% at higher masses starting at about 220 GeV for £ ranging from 0.002 to 0.1. In particular, this 
means that the low mass objects (Higgs bosons or top squarks) are mainly detected in the 420 m pots whereas the 
heavier ones in the closer pots at 210 m 6 . 

Our analysis does not assume any particular value for the £ resolution. We will discuss in Sections III B and IV C 
how the resolution on the Higgs boson or the top squark quark masses depend on the detector resolutions, or in other 
words, the missing mass resolution. 

III. SUSY HIGGS BOSON PRODUCTION 
A. Theoretical aspects 

Let us briefly recall the properties of the lightest MSSM Higgs boson h (for recent reviews, see Refs. [2^ and p39|). 
As is well-known, h is constrained to be lighter than the Z boson at tree-level. Once radiative corrections are taken 
into account 0, HE H^, the upper limit on its mass becomes mh ^135 GeV. The actual value of nih depends 
on several MSSM parameters: two parameters that are sufficient to describe the Higgs boson sector at tree-level, 
generally chosen to be itia and tan j3, the ratio of the vacuum expectation values of the two Higgs boson doublets of 
the MSSM; and additional parameters that control the size of the radiative corrections. These are the top squark and 
bottom squark soft supersymmetry breaking masses, assumed in this paper to be degenerate and denoted by Msusy, 
the top squark and bottom squark triscalar couplings (A-terms) A t and Ab, and the supersymmetric Higgs boson 
mass parameter \x. The dependence of the lightest Higgs boson mass on these parameters can be roughly described as 
follows: rrih increases with uia and tan/3, as well as with the common third generation squark mass Msusy- Its value 
also depends strongly on the top squark mixing parameter X t = A t — /icot/?: starting from the "minimal mixing" 
X t = 0, it increases with X t and reaches a maximum for X t ~ \^6Ms, where M| = (m| + m 2 )/2 is the average of 
the two top squark squared masses (Ms — Msusy in the limit Msusy S> rrit). This is illustrated by the following 



5 The expression of C h, C t t are explicitely given in formulae (2.13) and (2.14) of Ref.Qj and formulae (4) and (15) from Ref. Q, and 
included in DPEMC Bl . 

6 Information from the 220 m pots can be included into the first level of the trigger. However, the information from the pots located at 
420 m come too late to reach the first level trigger, and can only be used in the second level trigger. Therefore, one must use either 
asymmetric events (one tag at 220 m and another one at 420 m, where the trigger is based on the information from the 220 m pot) or 
directly a trigger from the main ATLAS and CMS detector. The latter option is challenging and under study in both collaborations. 



approximate formula for the one-loop upper bound on the lightest Higgs boson mass, valid in the decoupling limit 
m A > m z and for m t X t < Mj plE^ : 



m\ < mf cos 2 2/3 - 



8n 2 r 
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12M| 



(6) 



In the minimal mixing case, mj, can reach an upper limit of about 120 GeV for Msusy % 1 TeV, while it can reach 
about 135 GeV in the maximal mixing case |39|. 

The couplings of h can significantly depart from those of the SM Higgs boson. In particular, its tree-level couplings 
to down-type and up-type fermions (normalised to the SM Higgs boson couplings) are given by: 
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sin(/3 — a) + cot /3 cos(/3 — a) (f 
sin(/3 — a) — tan/3 cos(/3 — a) (f 



up-type fermion) 
down-type fermion) 



(7) 



where a is the angle that diagonalises the squared mass matrix of the CP-even Higgs boson and defines the physical 
CP-even states h and H. As for the couplings to the gauge bosons, hZZ and hWW \ they are suppressed by a 
factor sin(/3 — a) relative to their SM values. In the decoupling regime tua ^> mz [23|, Hj], in which A, H and 
are all much heavier than h, | cos(/3 — a)\ < 0(m 2 z /m A ) <C 1 and therefore the couplings of the lightest MSSM 
Higgs boson h approach those of the SM Higgs boson (in particular, ghff — 1)- On the contrary, if tua ~ mz 
(or more precisely tua < m™ ax , where m™ ax is the maximal value can reach for fixed values of the squark 
parameters), | cos(/3 — a)\ ~ 1 and therefore h has significantly different couplings from those of the SM Higgs 
boson. In particular, at large tan/3, in the "anti-decoupling" regime, its couplings to down-type fermions are strongly 
enhanced (|pft,6& | — \9hrr\ — tan/3 ^> 1), while its couplings to up-type fermions and gauge bosons are suppressed 
{\ghtt\ ~ cot/3 <C 1 and guww = ghzz = sin(/3 — a) <C 1, in units of the SM Higgs boson couplings). As we shall see 
below, this enhances the production cross-section of the lightest Higgs boson via gluon fusion, while the associated 
production with gauge bosons, qq — > Zh/Wh, is suppressed. Also the partial decay width of h into bb (t + t~), 
which is proportional to g\ hh (g 2 TT ), is enhanced. By contrast, the decay h — > 77, which in the decoupling regime is 
dominated by the W boson loop, does not benefit from such an enhancement at large tan (3 (the subdominant bottom 
quark loop is enhanced, but the dominant W boson loop is suppressed). This decay has therefore a suppressed 
branching ratio in the anti-decoupling regime. We close this short review of the anti-decoupling regime by noting 
that the heavier CP-even Higgs boson H, in contrast to h, has SM-like couplings, but is much heavier than h and A. 



Finally, another regime of interest is the "intense-coupling" regime |25j, which occurs when tua 



mi 



and tan /3 



is large. In this regime, all three neutral Higgs bosons are very close in mass, nih w tua ~ mn, and have enhanced 
(suppressed) couplings to down-type fermions (down-type fermions and gauge bosons - the couplings AWW and 
AZZ are forbidden by CP invariance), so that it may be difficult to distinguish among them at the LHC. 

Let us now discuss the production of the lightest MSSM Higgs boson via gluon fusion |i l li M IM ItR l4j | IZ9I I50I l5l) 
(see Appendix lAlfor the relevant formulae). In the SM, top quark loops give the main contribution to the cross-section, 
and bottom loops give a smaller contribution. In the MSSM, the contribution of the bottom loops can become very 
large at large tan/3 (in the regime where the hbb couplings, Eq. (|A5(I . are enhanced) while the top quark loops are 
suppressed, resulting in an enhancement of the gluon fusion cross-section. In addition, top and bottom squark loops 
contribute. However, top squark loops significantly affect the cross-section only in the case of a light top squark, 
TO ti S (200 — 400) GeV. In the decoupling regime, their effects are particularly spectacular in the presence of a large 
top squark mixing, in which case they interfere destructively with the top quark contribution p2|. For bottom squark 
loops to be sizable, a large value of tan/3 is also needed, as well as a large value of \/j,\ in the decoupling regime. 

In the regime we are interested in, which is characterised by a large value of tan/3 and a suppressed value of 
sin(/3 — a), there is no enhancement of the htt couplings at large top squark mixing, contrarily to the situation in the 
decoupling regime. However, the hbb couplings are enhanced at large mixing; but this is compensated by the fact that 
the bottom squark masses are generally larger than the lightest stop mass mi in typical MSSM scenarios. Therefore, 
we neglect the squark loops in the following discussion, although they are included in our numerical results. Neglecting 
as well the terms suppressed by cot /3, we then find the following enhancement factor for the MSSM cross-section with 
respect to the SM cross-section (the QCD corrections to the leading order cross-sections are expected to reduce this 
ratio |2y| by some 30% at large tan/3): 



VMssAijgg — > h) 
vsm (gg -> h) 



sin(/3 — a) — tan/3cos(/3 — a) 



A?(T t )+A%( n ) 



(8) 



where the loop functions A^(t) and A%(t) are defined in Eq. I|A2J| . and r t ^ = m\j '4mL. We therefore expect a 
large enhancement factor at large values of tan/3 in the regime where uia < m™ ax . This can indeed be seen in Fig. 



n (upper plot), where &MSSAi(gg — ► h)/asM(gg —> h) is shown as a function of rrih for tan/3 = 30 and various values 
of the squark parameters. In the maximal mixing case, m™ ax is large and the antidecoupling condition rriA < m™ ax 
is satisfied over the range 90GeV < rrih < 120 GeV (remember that in this regime rrih w tua), hence cos 2 (/3 — a) 
remains very close to 1, and the curve essentially reflects the dependence of the loop functions A^fa) and A^(tj,) 
on rrih- In the minimal mixing case, ni h aax is smaller, especially for Msusy = 500 GeV (namely rri!^ ax w 114 GeV 
for Msusy = 1 TeV, and m™ ax w 107 GeV for Msusy = 500 GeV), so that one leaves the anti-decoupling regime 
for much lower values of rrih than in the maximal mixing case. This explains why the enhancement factor strongly 
decreases when rrih approaches m™ ax , and finally reaches the decoupling regime value &MSSM(gg — ► h) / '<jsm{99 — * 
h) = 1 for rrih — rn™ ax (up to squark loop effects, which remain small for the squark parameters considered here). 
Fig. H (lower plot) shows the dependence of the enhancement factor on tan/3 for rrih — 100 GeV. For this value of 
rrih, the condition wia < rri h nax is satisfied for all four sets of the squark parameters considered. For moderate values 
of tan/3 (tan/3 < 5 is excluded experimentally for rrih — 100 GeV), the anti-decoupling regime is not yet reached, i.e., 
cos 2 (/3 — a) is large but not maximal (cos 2 (/3 — a) ~ 1). For larger values of tan /3, e.g., tan/3 > 20, the anti-decoupling 
regime is reached and the enhancement factor grows with tan 2 /3, as expected. 

B. SUSY Higgs boson production at the LHC 

In this section, we address the MSSM Higgs boson production for masses below 120 GeV, when the Higgs boson 
decays into bb, the least favourable case at the LHC. As mentioned in the previous section, Figure ^ shows the cross 
section enhancement factor for SUSY Higgs boson production with respect to the Standard Model case at generator 
level. In the upper plot of Fig. ^ the full and dashed lines show the results for the minimal mixing scenario for two 
common values of the mass, Msusy, of third generation MSSM squarks (1000 and 500 GeV). They lead to typical 
masses of the top squark and bottom squarks of 1010 or 520 GeV, respectively. The cross-section was computed using 
bottom, top, top squark and bottom squark loops, while the effect of top squark and bottom squark loops is less than 
one per mil. The enhancement factor can go up to a factor 20 compared to the Standard Model case, but is very 
dependent on the mass of the Higgs boson. In the maximal smearing scenario (dotted and dashed dotted curves in 
Fig. Q, the enhancement factor is found to be similar to that at low Higgs boson masses, but remains important at 
higher masses. 

The bottom plot of Fig. displays the dependence as a function of tan/3 for a Higgs boson mass of 100 GeV and 
for the same scenarios as before. The enhancement factor for the Higgs boson production cross section can reach a 
factor of up to 45 for a value of tan/3 of 50. For this particular value of the mass of the Higgs boson, the model 
dependence is not very large. 

It is important to note that searches will benefit directly from the increase of the cross-section since they will be 
looking for Higgs bosons decaying into bb in the main detector, and the branching ratio of h — ► bb is quite stable as 
a function of the MSSM parameters in this region of phase space 7 . Thus, the search for diffractive production of 
MSSM Higgs bosons is the only one benefiting fully from the increase of the cross-section at high values of tan /3. 

In the following, we perform a detailed study of signal-over-background ratio in the case of a Higgs boson mass of 
120 GeV. We chose this particular mass since in most of the MSSM parameter space, the Higgs boson mass is below 
this value, and this mass leads to the least favourable scenario (the lowest cross-section and signal-over-background 
ratios) with respect to lower masses. In Fig. we give the signal-over-background ratio for Standard Model and 
MSSM Higgs boson production for a mass of the Higgs boson of 120 GeV and for different values of tan /3, as a function 
of the roman pot mass resolution (corresponding to the Higgs boson mass resolution) for a luminosity of 100 fb _1 . 
This study was performed after a fast simulation of the CMS detector (the ATLAS detector simulation is expected 
to produce very similar results) and experimental cuts described in the following paragraph. 

First of all, we require both final state protons to be detected in the roman pot detectors, and we take into account 
the acceptance of these detectors as it is discussed in section I. The cuts applied in the analysis are detailed in Ref. 0. 
The basic idea is to require two high pr b-jets with pn > 45 GeV, pti > 30 GeV, originating from the decay of the 
Higgs boson into bb at low masses. The difference in azimuth between the two jets should be 170 < A$ < 190 degrees, 
asking the jets to be back-to-back. Both jets are required to be central, \rj\ < 2.5, b-tagged, with the difference in 
rapidity of both jets satisfying |A^| < 0.8. A cut is applied on the ratio of the dijet mass to the total mass of all 
jets measured in the calorimeters, Mjj/M a u > 0.75. The ratio of the dijet mass to the missing mass should fulfil 
Mjj/faba) 1 ' 2 > 0.8. 



7 This is not the case when one looks into non diffractive MSSM Higgs bosons decaying into 77 which is strongly suppressed at high 
tan/3, see paragraph III A. 




FIG. 1: Enhancement factor for the diffractive Higgs boson production cross section. Upper plot: enhancement factor as 
a function of the Higgs boson mass for a value of tan/3 = 30, and different mixing scenarii and SUSY masses (full line: 
minimal mixing, Msusy = 1000 GeV; dashed line: minimal mixing, Msusy = 500 GeV; dotted line: maximal mixing, 
Msusy = 1000 GeV; dashed dotted line: maximal mixing, Msusy = 500 GeV). Lower plot: similar study as a function of 
tan (3 for a Higgs boson mass of 100 GeV. The MSSM Higgs boson spectrum has been obtained using the program FeynHiggs 
with n = 200 GeV, M 3 = 800 GeV, M 2 = 200 GeV and Mi = (a 1 /a 2 )M 2 . 



The case for the Standard Model Higgs boson was already given in |3j , and we follow the same approach concerning 
the background and signal studies. To compute the signal over background ratios, both signals and backgrounds 
dominated by exclusive bb production have been integrated over a 2 GeV mass window. After cuts, the typical 
number of events expected for the signal of a 120 GeV Higgs boson and for a luminosity of 100 fb _1 is 27.1, 73.2, 154, 
398 and 1198 for Standard Model and MSSM (tan/3= 15, 20, 30 and 50) Higgs boson production. If the Higgs boson 
is supersymmetric and if tan [3 is large, the diffractive production of MSSM Higgs bosons could lead to a discovery in 
the double pomeron exchange mode at the LHC. Figure demonstrates that the signal over background can reach a 
value up to 54, 26, 16, and 13 for respective Higgs boson mass resolutions in roman pot detectors of 1, 2, 3 and 4 GeV 
and for a value of tan/3 of 50 for a luminosity of 100 fb _1 . 
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FIG. 2: Signal over background as a function of roman pot mass resolution for a Higgs boson mass of 120 GeV and for different 
values of tan/3. From bottom to top: full line: Standard Model Higgs boson, dashed line: SUSY Higgs boson with tan/3 =15, 
dotted line: tan/3 =20, dashed dotted line: tan/3 =30, full line: tan/3 =50. 



IV. PRODUCTION OF TOP (BOTTOM) SQUARK PAIRS 
A. Theoretical framework 



In the MSSM, for each quark flavour q, there are two supersymmetric scalar partners qx and qn associated with 
the two fermion chiralities q^ and qR. In general, these scalars are not mass eigenstates, due to the presence of soft 
supersymmetry breaking terms which mix them, the A-terms A q y q QLq R H + h.c, where y q is the Yukawa coupling 
of the quark q and H is one of the two MSSM Higgs boson doublet. The mass eigenstates q\ and q2 are obtained by 
diagonalising the following 2x2 matrix 55] : 



M 2 Q 



Dr 



m q X q 



Ml 



m q X q 



Dj 



(9) 



where Dl = (T q — Q q sin 2 Qw)m\ cos 2/3, Dr = Q q sin 2 6wm% cos 2/3, X q = A q — ficotfl for up-type squarks and 
X q = A q — /itan/3 for down-type squarks. The soft-supersymmetry-breaking squark masses Mq and Mr are of the 
order of the supersymmetry-breaking scale Msusy, and phenomenological constraints require the A-tevm parameter 
A q to be at most a few times Msusy |56j |. Neglecting the terms proportional to m 2 z in the squark mass matrix, the 
mass eigenvalues and the mixing angle, which relates the weak interaction eigenstates qu,R to the mass eigenstates 
91,2, are given by the following expressions: 



91,2 



= -( 



Mi 



m T 



(Mq - M R ) 2 



tan Op; = 



2m q X q 



Mi 



Ml 



(10) 



In practise, the mixing is significant only if m q X q ~ Mg USY , which can occur for the top squark and, at large tan/3, 
also for the bottom squark. In this case, one can have a strong hierarchy between the two mass eigenstates, rriq 1 -c m ? - 2 . 
The 95% C.L. experimental bounds on the lightest top squark and bottom squark masses are mj > 95.7 GeV and 
mg > 89 GeV, respectively, while the bound on the other squarks is 250 GeV |5j| (the latter bound also applies 

to b\ if mixing effects are small in the bottom squark sector). Although these bounds where derived under specific 
assumptions and may therefore not hold in some regions of the MSSM parameter space, they are rather robust. 



In our experimental study, we consider the following values for the lightest top squark mass: 174.3 GeV (i.e., 
= mt), 210 GeV and 393 GeV. As we will see, the resolution that can be obtained on the top squark mass 
crucially depends on its decay width, which in turn is a function of the top squark mass and of the other MSSM 
parameters. If the top squark is very light, it is likely to be the next-to-lightest supersymmetric particle, assuming the 
lightest supersymmetric particle (LSP) is the lightest neutralino Xi- Then all two-body decay channels occurring at 
the tree level are closed. If in addition the tree-body decay channels t\ — > bW + x.i and t\ — > bH + Xi are kinematically 
not accessible, the main decay mode is expected to be the loop-induced flavour violating decay i\ — > cxi .58]. The 
decay width of the lightest top squark is then given by [H^ : 

„2 / m 2 \ 2 

r(*i-cx?) = ^l/£xl 2 N 2 m fl ^-^fj . ( n ) 

where \ff J - i \ < 1 is the Cl-Cl-Xi coupling, and e is a flavour-violating insertion. The authors of Ref. (5^ estimated 
|e| ~ (1 - 4) X 10~ 4 in mSUGRA, yielding T(ti -► ex?) < (0.085 - 1.4) x 10~ 9 m h [1 - (m~+ /m^) 2 } 2 ; but depending 
on the mSUGRA parameters |e| could be either much smaller or larger, in particular at large tan/3 where |e| behaves 
like tan 2 (3. In non-minimal SUGRA models, |e| could even be of order one. However, in the regions of the MSSM 
parameter space where T(t\ — > cxi) is suppressed, the four-body decay modes t\ — > bxiff are likely to be dominant 
For larger top squark masses, the three-body decay channels i\ — > bW + Xi and t\ — > bH + Xi |6lll62j (and, if 
the sleptons are lighter than the lightest top squark, t\ — > bv\ l + and i\ —> bi>i l + jEEE3)E3) °P en U P and tend to 
dominate. Finally, when > mb + and > m t + m^o , the two-body decays t\ — » bx^ and t\ — > txl become 

kinematically accessible and dominate the lightest top squark decays 8 (64|. The partial decay width of ?i — > bx\ is 
given by: 

r(*i -> bxt) 

where In and k\\ are chargino couplings. Since \k\\\ < 1, T(b\ — * bxt) ~ 0.0085 mf t [l — (m-^/m^) 2 ] 2 . The 

decay t\ — > tx\ has a larger phase space suppression. 

For our experimental study, we do not consider any specific benchmark scenario, but simply assume < 
for the cases m fl = 174 and 210 GeV. Then, we conservatively take T it = 100 MeV for m~ ti = 174 and 210 GeV, 
although the actual decay width could be much smaller. For = 393 GeV, we assume that the two-body decay 
t\ — > 6xi" is accessible, and wc take the decay width computed for SPS la, If =1.8 GeV |29| . 



16tt 



(Z 2 1 +fc? 1 )m fl 1- 



2 

*i 



(12) 



B. Stop production cross-section and missing mass distribution 

At hadron colliders, top squarks can be produced at lowest QCD order via quark-antiquark annihilation and gluon- 
gluon fusion. In the present study, we are interested in the J z = 0, colour-singlet gluon-gluon fusion cross-section; it 
reads, at the parton level plj : 

d _ZML {gg ^~ ti k) = (13) 

dt Vyy ' 12 s 2 K ' 

where ^fs is the centre-of-mass energy of the subprocess, m^. is the top squark mass, and Et is the transverse energy 
of the final particles. 

The top squark production cross-section has been obtained using the DPEMC generator [22] after applying a survival 
probability of 0.03. The top squark pair production cross-section as a function of the top squark mass is given in Fig. 
13 The tt production cross section is found to be 26.3, 14.1 and 1.1 fb for a top squark mass of 174.3 (at about the 
top quark mass), 210 and 393 GeV respectively. 

The distribution of the missing mass distribution for tt and tt events for = m t =174.3 GeV is shown in Fig. 0] 
for a luminosity of 100 fb" 1 . The missing mass distribution for top squark (top) events is in full (dashed) line. The 



We do not consider the decays t\ — > tg and t\ — > biW + , biH + , since the gluinos and the other squarks are generally heavier than the 
lightest top squark. 



cross-section rise at threshold is much faster than for top quarks and typical of pair production of scalar particles. 
The next section describes a method to determine the stop quark mass by performing a threshold scan of the missing 
mass of the tt process, measured with roman pot detectors. 

C. Stop mass measurement 

In this section, we describe briefly the method we used to obtain the stop mass resolution and its results. The 
histogram method 9 is described in more detail in Ref. |37l |. It compares the mass distribution in data with some 
reference distributions following a Monte Carlo simulation of the detector with different input masses corresponding 
to the data luminosity. As an example, we can produce a data sample for 100 fb _1 with a top squark mass of 210 GeV, 
and a few MC samples corresponding to top squark masses between 180 and 240 GeV by steps of 1 GeV. To evaluate 
the statistical uncertainty due to the method itself, we perform the fits with some 100 different "data" ensembles. For 
each ensemble, one obtains a different reconstructed top squark mass, the dispersion corresponding only to statistical 
effects. The \ 2 1S defined using the approximation of poissonian errors as given in Ref. j3^- Each ensemble thus gives 
a x 2 curve which in the region of the minimum is fitted with a fourth-order polynomial. The position of the minimum 
of the polynomial gives the best value of the top squark mass and the uncertainty cr(mi) is obtained from the values 
where x 2 = Xmi„ + 1- 

The results are given in Fig. [5] and© Fig. ©displays the results on the top squark mass resolution for a top squark 
mass of 210 GeV as an example, as a function of luminosity, for different roman pot resolutions. The results depend 
only weakly on the roman pot resolution and mostly on the number of events produced for a given luminosity. The 
resolution on the top squark mass is thus dominated by statistics. We also note that the integrated luminosity does 
not take into account the efficiency of the cuts to select the tt events since these efficiencies depend strongly on the 
SUSY parameters. A typical efficiency of 60% is found requesting a missing transverse energy to be greater than 80 
GeV, and either two reconstructed jets or one lepton and one jet with a transverse momentum greater than 20 GeV. 
Figure © displays the resolution obtained for the three values of the top squark mass discussed above. A resolution of 
about 0.4, 0.7 and 4.3 GeV is obtained for a top squark mass of 174.3, 210 and 393 GeV for a luminosity (divided by 
the signal efficiency) of 100 fb" 1 . As it was mentionned in paragraph IV A, the top squark width has been taken into 
account in this study. For a top squark mass of 174.3, 210 GeV, the top squark width of 100 MeV has a negligible 
effect, whereas the top squark width of 1.8 GeV for a top squark mass of 393 GeV cannot be neglected. 

V. CONCLUSION AND OUTLOOK 

In this paper, we described the advantages of diffractive SUSY particle productions for two different processes, 
namely the MSSM Higgs bosons and top squark pairs. The cross-section for diffractive MSSM Higgs bosons production 
is noticeably enhanced at high values of tan (3 and since we look for Higgs boson decaying into bb, it is possible to benefit 
directly from the enhancement of the cross section contrary to the non diffractive case. A signal-over-background up 
to a factor 50 can be reached for 100 fb _1 for tan/3 ~ 50. In particular, we analysed in detail the antidecoupling 
regime, in which H behaves like the SM Higgs boson, while h has enhanced (resp. suppressed) couplings to down- 
type fermions (resp. up- type fermions and gauge bosons). We find that central diffraction production seems to be 
promising in that regime. 

The other application is to use the so-called "threshold-scan method" to measure the top squark mass in exclusive 
events. The idea is quite simple: one measures the turn-on point in the missing mass distribution above twice the top 
squark mass. After taking into account the top squark width, we obtain a resolution on the top squark mass of 0.4, 
0.7 and 4.3 GeV for a top squark mass of 174.3, 210 and 393 GeV for a luminosity (divided by the signal efficiency) 
of 100 fb^ 1 and the production rates calculated in Section IV of this paper. If these rates hold, the typical mass 
resolutions are comparable to those at a linear collider. The process is thus similar to those at linear colliders (all 
final states are detected), but without the initial state radiation problem. 

It should be stressed once more that production via the diffractive exclusive processes is model dependent, 
and definitely needs the Tevatron data to test the models. It will allow to determine more precisely the production 
cross-section by testing and measuring at the Tevatron the jet and photon production for high masses and high dijet 
or diphoton mass fraction. If, for instance, we compare (as much as possible since the methods of evaluation are not 



In Ref. |37l . we give two methods to measure the W boson or the top mass, namely the histogram or the turn-on fit methods. For a 
matter of simplicity, we used only the histogram method in this paper. 
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FIG. 3: Top squark pair production cross-section as a function of the top squark mass. 
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FIG. 4: Distribution of the missing mass for 100 fb _1 for ti events (dashed line), and for it (full line) for mt op = m^. The faster 
rise of the stop quark cross-section as a function of missing mass is due to the scalar nature of these particles. 



exactly comparable) the expectations of the 'Pomeron-induced' 0, Q and 'proton-induced' Q models for exclusive 
production, we find striking differences for high mass states, in particular the stop quark production. Indeed, the 
Sudakov form factors present in the 'proton-induced' approach induce negligible stop quark cross-sections l21| . The 
search for exclusive events in high mass dijet or diphoton events at the LHC is of considerable importance to test the 
models. For MSSM Higgs bosons in the low mass range there is not an order-of-magnitude difference between the 
two models as was already mentioned for the standard Higgs production [lfij . 

On the other hand, it is also possible to perform a similar study using inclusive double pomeron exchanges (J2j. 
These processes have already been observed by many experiments but suffer from the lack of knowledge on the gluon 
density in the pomeron at high (3. The first step is thus to measure this gluon density by, for instance, using dijet 




FIG. 5: Expected statistical precision of the t mass as a function of the integrated luminosity for various resolutions of the 
roman pot detectors using the histogram-fitting method. 
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FIG. 6: Expected statistical precision of the t mass as a function of the integrated luminosity for different t masses (174.3, 210 
and 393 GeV). 



events or the threshold scan method for inclusive ti production. Once the high-/? gluon better determined, it is 
possible to look for top squark events, again using the threshold scan method and deviation at high masses provided 
the cross-section is hig h enough. This study goes beyond the purpose of the present paper and it certainly deserves 
a dedicated study |37| . 
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APPENDIX A: MSSM HIGGS BOSON PRODUCTION VIA GLUON FUSION AT LEADING ORDER 



The cross-section for the lightest MSSM Higgs boson production via gluon fusion reads, at leading order |5j 

2 



a(gg -> h) 



Gfo^ 
288V2tt 



4 '— J m~ 

q 1 



(Al) 



where the first term contains the quark loop contributions, and the second term the squark loop contributions. The 
loop functions A^(t) and A^(t) are given by: 



A h g (r) 



1[t + (t-1)/(t)], 4(r) 



[f(r)-r], 



(A2) 



arcsin 2 ( v /r) 



fir) 



ln ,i±V^ 



T < 1 

r > 1 



(A3) 



and the parameters r q and Tq are defined by r q = m 2 /4m 2 and = m^/4m~, respectively, where m q (resp. m q ) 
denotes the mass of the quarks (resp. squarks) running in the loop. The couplings of the lightest MSSM Higgs boson 
to the top and bottom quarks, normalised to the SM Higgs boson couplings, are given by: 



ghtt = sin(/3 — a) + cot (3 cos(/3 — a) , 
9hbb = sin(/3 — a) — tan (3 cos(/3 — a) , 



(A4) 
(A5) 



and its couplings to the top squark and bottom squark mass eigenstates, in units of g/Myy, by (we omit the off-diagonal 
couplings g^i^ and g h i i , which are not relevant at leading order): 
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(A6) 
(A7) 
(A8) 
(A9) 



where 9^ and 9~ h are the mixing angle in the stop and the bottom squark sector respectively, defined by q\ = cos 9 q ql + 
sin 9 q q R , q 2 = - sin 9 q q L + cos 9 q q R . 

The leading order cross-section for the SM Higgs boson production via gluon fusion can be obtained from Eq. (|A1|) 
by removing the squark contribution and by setting gh qq = 1. 
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